Cloning quantum entanglement in arbitrary dimensions 
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We have found a quantum cloning machine that optimally duplicates the entanglement of a pair 
of d-dimensional quantum systems. It maximizes the entanglement of formation contained in the 
two copies of any maximally-entangled input state, while preserving the separability of unentangled 
input states. Moreover, it cannot increase the entanglement of formation of all isotropic states. For 
large d, the entanglement of formation of each clone tends to one half the entanglement of the input 
state, which corresponds to a classical behavior. Finally, we investigate a local entanglement cloner, 
which yields entangled clones with one fourth the input entanglement in the large-d limit. 



I. INTRODUCTION 

The no-cloning theorem Q] precludes a perfect cloning 
of arbitrary quantum states. However, an imperfect 
cloning is possible and various quantum cloning machines 
(QCM), which duplicate quantum states with the high- 
cst fidelity, have been proposed following the seminal pa- 
per of Buzek and Hillery Q. Recently, the question of 
whether quantum entanglement can be cloned or not was 
raised in Q. Since the quantum entanglement is a re- 
source for quantum computation, quantum communica- 
tion, and quantum cryptography, it is important to know 
up to what extent this resource can be duplicated. For 
maximally-entangled (ME) states in two dimensions, an 
entanglement no-cloning principle was formulated : "if 
a quantum operation can be found that perfectly du- 
plicates the entanglement of all ME states, then it is 
necessary does not preserve separability". A QCM was 
proposed that optimally (but imperfectly) clones the en- 
tanglement of two-dimensional bipartite states (qubits) 
while preserving separability. 

In the present paper, we extend these results to pairs 
of (i-dimensional quantum systems, with arbitrary d. We 
show that a (symmetric) cloning machine can be defined, 
which maximizes the amount of entanglement of the two 
clones of ME states, while producing separable clones in 
the case of unentangled input states. We analyze the en- 
tanglement of the clones in terms of fidelity, but show 
that optimizing the cloning machine in terms of fidelity 
actually leads to maximizing the entanglement of forma- 
tion of the clones provided that we restrict to cloning ma- 
chines that are covariant under local unitàries. We then 
compare the resulting optimal d x d entanglement cloner 
to a "local" cloning transformation that can be achieved 
by applying a separate universal cloning machine to each 
component of the bipartite system. 



and (ii) maximizes the entanglement of the two clones 
resulting from any ME input state. We will characterize 
a cloner by considering the transformation of an input 
that is maximally entangled with a reference system. By 
projecting the reference onto (the complex conjugate of) 
the input state, one gets the operation of the cloner on 
this state. The general form for such a cloning transfor- 
mation is defined in the computational basis {\i}} by the 
state 



\S)n,a,b,A= ^2 s i]ki\i)n\j)a\k)b\l). 



(1) 



where 1Z denotes the reference system, a and ò stand for 
the two clones, and A corresponds to the ancilla. All the 
summations here are d 2 -dimensional since all the states 
involved are <i 2 -dimensional bipartite states, e.g., \i) — 
\ía)\Íb}- Thus, each index i, j, fc, or l actually represents 
a couple of indices running each from to d — 1, e.g., 
i — {ía, with ía-, %b € [0; d — 1] . Of course, the index 
A stands for Alice's component of the bipartite states, 
while B stands for Bob's component. 

As mentioned above, the joint state of the two clones 
and the ancilla is obtained by performing an appropriate 
projection on the reference system. Thus, for an input 
state |<£>) = X)i n i|*)i the result of the cloning transfor- 
mation is of the form 



i,j,k,l 



II. COVARIANT CLONER UNDER LOCAL 
UNITÀRIES IN DIMENSION d x d 

Following the ideas presented in Q, we seek for a 
cloning transformation that (i) preserves separability, 



Then, the state of any one of the clones is further ob- 
tained by tracing out the ancilla and the other clone. 
This entanglement cloning transformation is shown in 
Fig. 1. 
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FIG. 1: (i) Optimal dxd entanglement cloner; (ii) "local" entanglement cloner, as defined in Sec. IV.Here, A and B stand for 
Alice's and Bob's part of the bipartite state, a and b stand for two clones. Entanglement is indicated by double loops. 



Next, we impose the following covariance condition on 
our cloning machine. Since we know that any local uni- 
tary operation acting on the A and B components of 
a bipartite state preserves its entanglement, we require 
that any such transformation acts similar ly on the clones. 
This condition amounts to imposing 

\S)kM = U*®U®U®U* \S) n ,a,b,A (3) 

where U is the product of any two unitary transforma- 
tions acting separately on each c?-dimensional component 
of the bipartite state, that is 

U = U A ® U B (4) 

where the indices A and B denote Alice's and Bob's com- 
ponents. Defined in this way, the operator U possesses 
a SU(d) (g> SU(d) symmetry. The covariance property 
implies that Sijki in Q is an invariant tensor of rank 4, 
which satisfies 

s ljk i = U^Ujj'Ukk'U^sefk'i' (5) 

where U* denotes the matrix element-wise complex con- 
jugate of U with respect to the computational basis. 
(Here, the summation over all repeated indices is im- 
plícit.) The matrix elements £/y = (i\U\j) satisfy 
U* k Ukj — Sij. If the transformation is real (i.e. Uij — 
U*j), then the invariant tensors for Alice's component 
are of the general form 5 tAjA ó kA i A , or S ÍA i A S jA k A , or 
5i A k A Sj A i A (see Q). For a complex transformation on 
Alice's component, the terms of the form Si A i A Sj A k A are 
excluded, so we are left with the two remaining terms. 
Then, for a tensor product of (complex) unitary trans- 
formations on Alice's and Bob's components, a general 
form for the invariant tensor is 

Sijkl 

= As i A 3A^k A i A 5 lB]B 5k B i B + B S ÍA j A Sk A i A S ÍB k B Sj B i B 
+ C ' S ÍA k A Sj A i A S ÍB k B Sj B i B + D 8 ÍAkA 5 jA i A 5 ÍBjB 8k B i B . 

(6) 



For a symmetric cloner, the output state must be in- 
variant under the interchange of the two clones, i.e., un- 
der permutations {íatÍb) *-* Q^a^b)- This implies that 
A = C and B = D, so we are left with only two param- 
eters A and B to be determined. 



III. OPTIMAL ENTANGLEMENT CLONER IN 
DIMENSION dxd 

The covariance condition J3J guarantees that our QCM 
transforms all states which are equivalent up to local uni- 
tàries (which have therefore the same entanglement) into 
equally entangled clones. In particular, the clones of all 
ME states will be equally entangled. Then, a cloner that 
is optimized on a particular ME input state will be opti- 
mal for all ME states. We choose as an initial dxd ME 
state 

l*)= J2 n ÍAÍ B \ÍA)\ÍB) (7) 
i A ,i B =0 

where rii A .i B = &í a ï b /y/d. As we shall show later on, we 
can maximize the entanglement of the clones simply by 
optimizing our QCM in terms of the fidelity of the clones, 

F = <#|p a |$), (8) 

where 

Pa = Tr A>b [\ X )(x\] (9) 

is the state of clone a. For the ME state 0, this fidelity 
is found to be 

F = \A\ 2 (d 2 + 3) + 4|S| 2 + 4ÏÏ(AB*)^-^. (10) 

Taking into account the normalization condition for the 
joint output state |%), 

2 (|^1| 2 + \B\ 2 ) (d 2 + 1) + 8dïï(AB*) = 1 , (11) 
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we can maximize the fidelity, Eq. (|10fl . which yields 
I / d 2 + 1 



F = 




4 \ d 2 - 1 



<P_T 
d 2 \d 2 -l 



(12) 



Note that, for d = 2, this result coincides with the max- 
imal fidelity of the entanglement cloner for two qubits 
obtained in 3], namely 



F 



12 



0.7171. 



The corresponding solution in d dimensions is 
.4 



dy/1 + Y(d) -y/1- Y(d) 



2(d 2 - 1) 



B = — 

where 

Y(d) = (l- 



dy/1 - Y{d) + y/1 + Y{d) 

2(d 2 - 1) 
(d 2 - 2f 



d 2 (d 2 - l) 2 +4(d 2 - 2) 2 



1/2 



(13) 

(14) 
(15) 

(16) 



IV. COMPARISON WITH OTHER CLONERS 

We compare the fidelity achieved by our optimal dxd- 
dimensional entanglement cloner, Eq. 112JI . with that of 
the universal cloner p| 



F u - 2 



1 



d 2 + l'- 
as well as that of the optimal real cloner 



F T 



Vd 4 + 4d 2 + 20 -d 2 + 2 
4(d 2 + 2) ' 



(17) 



(18) 



In order to make this comparison consistent, we have ob- 
tained formulac (|17|l and Ijl8(l by replacing the argument 
d by d 2 in the original formulae. This is done because, 
in our consideration, the dimension d stands for the di- 
mension of each component of the bipartite input state, 
so that the total dimension of our input state is d 2 . 

In Table I, we compare the fidelity F of our entan- 
glement cloner with F u and F r , for several vàlues of the 
dimension d. Our cloner performs better than the univer- 
sal cloner in d 2 dimensions for all d, which is obviously 
due to the fact that the ME states span only a subset 
of the entire set of cP-dimensional states. In contrast, 
the real d 2 -dimensional cloners outpcrform our cloners, 
exccpt if d = 2 where they coincide [3j. This can be in- 
terpreted by noting that the set of c? 2 -dimensional real 
states is generated by SO(d 2 ), with (d 2 (d 2 + l)/2) - 1 
real degrees of freedom, while the set of ME states is 
generated by SU(d) x SU(d), with (d 2 - l) 2 real degrees 
of freedom. For d — 2, they coincide, so that our cloner 



d x d 


F r 


F 


F u 


Floc 


2x2 


0.7171 


0.7171 


0.7000 


0.5833 


3x3 


0.6069 


0.6019 


0.6000 


0.4583 


4x4 


0.5617 


0.5592 


0.5588 


0.4000 


5x5 


0.5385 


0.5386 


0.5385 


0.3667 


6x6 


0.5277 


0.5271 


0.5270 


0.3452 



TABLE I: Optimal fidelity F of the dx d entanglement cloner 
for various dimensions d. It is compared to the fidelity of the 
real cloner F r and universal cloner F u , both in d 2 dimensions, 
and to the fidelity of the dxd "local" cloner Fi oc . The fidelities 
are shown in decreasing order. 



provides the same fidelity as that of the real cloner in 
dimension 4, namely Eq. (|13|l . This is related to the fact 
that the set of ME 2-qubit states is isomorphic to the set 
of 4-dimensional real states 0- For d > 2, the set of ME 
states is in some sense "larger" than the set of real states 
for d > 2, so that the achievable fidelity of the entan- 
glement cloner is lower. The fidelity of our cloner drops 
faster than that of the real cloner with increasing d, but 
remains always higher than the fidelity of the universal 
cloner. As expected, in the limit d — > oo, all three fideli- 
ties tend to the asymptotic value 1/2. In this limit, all 
quantum cloners can be interpreted simply as a classical 
transformation that maps the original state to one of the 
clones, chosen at random, the other clone being prepared 
in a maximally mixed state. 

Interestingly, we may also compare Eq. (|12fl to the fi- 
delity of a "local cloner" obtained by applying a cloner 
separately to Alice's and Bob's components (see Fig. 1). 
Since the state of Alice's or Bob's subsystem is maximally 
mixed (hence non-polarized) when the bipartite state is 
ME, it is natural to use a universal ci-dimensional cloner. 
We may observe that if we consider a cloning transfor- 
mation that performs such a local universal cloning, then 
it is represented by a joint state of the same type as Eq. 
0j see p. The only difference is that in the expression 
© for the tensor Sijki, all coefficients must be equal, i.e., 
A = B = C = D. The normalization condition (|ll·|) then 
gives immediately A = l/(2(d + 1)). Substituting this 
expression into Eq. I|10|l results in the fidelity 



d + 2 
2d{d+l) 



(19) 



for the local cloner. This fidelity is compared in Table 1 
with that of the other cloners. It appears that cloning 
Alice's and Bob's parts locally leads to a much lower 
fidelity. Note that for d — 2, the value of the fidelity 
Fioc in Table 1 coincides with the value 7/12 obtained in 
PJ. In the limit d — > oo, this fidelity tends to 1/4, which 
can be easily understood in the classical picture above. 
To contribute to the fidelity, both cloners indeed need to 
map Alice's and Bob's components of the original state 
onto the right clone, which only happens with probability 
(l/2) 2 = l/4. 
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V. ENTANGLEMENT OF FORMATION OF 
THE CLONES 

In order to investigate the entanglement properties of 
our cloning transformation, we will use as an entangle- 
ment measure for the clones the entanglement of forma- 
tion , which was computed for several classes of states 
that are invariant under some groups of local symme- 
tries In particular, we will be interested in the class 
of states that are invariant under the transformations 
U x U* for all U, called isotropic states in ja B - These 
states may be written in a general form as |10| 

p =L·J-(t-\<ï>)(<ï>\)+F\$){$\ (20) 

where 1 is the identity and |$) is given by Eq. J7J. Due to 
the covariance condition |J3J, our QCM transforms UxU* 
invariant states into into states that are also invariant un- 
der UxU*. We can check that, by cloning the particular 
ME state |<E>), which is U x U* invariant, we obtain a 
clone of the form 



Pa 



(d 2 + 2)\A{ 2 + 2\B\ 2 + Ad^(A*B)) |$)($ 



(\A\ 2 



2\B\ 



-S(A*5) 1 
d 



(21) 



which is indeed an isotropic state and is consistent with 
Eq. HI (jp . Hence, as a consequence of our covariance con- 
dition, all ME states, which can be obtained from | $) by 
applying local unitàries, are cloned onto isotropic states. 
For this class of states, the entanglement of formation is 
known for arbitrary dimensions [g, lïo| 



E F {p) 



0. 



■1 4(rf-l)i 



(F-l) + log a d, Fe[ 



~3T~ 



1] 

(22) 



where 

Rl, d -!(F) = ff 2 ( 7 (F)) + [l- 7 (F)]log 2 (d-l),(23) 



7 (F) = \ VF+ 



(24) 



H 2 (p) = -plog 2 (p)-(l-p)log 2 (l-p). (25) 

As shown in Fig. 2, the entanglement of formation Ep(p) 
is a monotonically increasing function of the fidelity F for 
isotropic states in any dimension d. Therefore, by opti- 
mizing our QCM in terms of fidelity we maximize, at the 
same time, the entanglement of the clones measured by 
their entanglement of formation. The circles in Fig. 2 
correspond to the maximal fidelity F that is achieved by 
our entanglement cloner, Eq. I|12f>. They show as well the 
corresponding entanglement of formation in this dimen- 
sion. The crosses mark the crossover between the expres- 
sion of the fidelity corresponding to the second and third 



lines of Eq. I|22|) . One notes that only for d > 7 there 
are vàlues of the fidelity for which the entanglement of 
formation has to be evaluated according to the third line 
of Eq. 1(22)1 . 




FIG. 2: Entanglement of formation Ef of the clone of a 
maximally-entangled input state versus the fidelity F of the 
clone for various dimensions d — 2 — 20 (the lowest curvc 
corresponds to d = 2, while the highest curve corresponds 
to d — 20). The circles show the maximum achievable fi- 
delity and the corresponding entanglement of formation. The 
crosses mark the crossover between the expression of the fi- 
delity corresponding to the second and third lines of Eq. 1221 1. 



In order to visualize how the entanglement itself is 
cloned, we plot in Fig. 3 the entanglement of formation 
Ef of the clones as a function of the entanglement of 
formation of the input ME state E- in , which is simply 
the von Neumann entropy of the reduced density ma- 
trix Ei n — log 2 d. We note that the entanglement of the 
clones is always less than one half the entanglement of 
the input state, while it asymptotically approaches this 
value for large d. The apparent "discontinuity" (if one 
can say so for a discrete graph) corresponds to d = 7, 
i.e., the crossover between the second and the third line 
of Eq. i|22|) when calculating the entropy of formation. 
In the limit of d — ► oo, the third line of Eq. (|22|l tends to 
Ef = F log 2 d = F E- m . Since the cloner can be viewed, 
in this limit, as a classical random distribution process 
associated with a fidelity F — 1/2, then the entangle- 
ment of each clone tends to one half the entanglement of 
the initial state. 

In Fig. 3, we also plot the entanglement of formation 
resulting from the "local" cloner discussed above. Re- 
cali that this cloner differs from our optimal (non-local) 
entanglement cloner only by setting A = B. Therefore 
it is also covariant, satisfying Eq. ©, and all our argu- 
ments about the entanglement of formation of the clones 
are applicable to this cloner as well. Thus, using the fi- 
delity of the clones (|19|) , we may plot the entanglement 
of formation of the clones. We see that the local cloner 
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FIG. 3: Entanglement of formation Ef of the clones of 
a maximally-entangled state obtained by the optimal (non- 
local) cloner (o) and the local cloner (+) versus the entangle- 
ment of the input state Ei n for various dimensions d — 2—200. 
The apparent "discontinuity" in both curves is due to the 
crossover from the second to the third line of Eq. 1221 for 
d > 7 (optimal cloner) and d > 13 (local cloner). Sòlid lines 
represent the asymptotics of Ef for large d in both cases. 



leads to a lower entanglement of formation, and even 
the asymptotics of Ep for large d is no more than one 
fourth the entanglement of the input state. The reason 
is that in the limit of large d, the classical random distri- 
bution only succeeds with probability 1/2 independently 
for Alice's and Bob's components, so the fidelity is 1/4. 
Hence, Ep — > E- m /A. These observations confirm that by 
increasing the dimensionality, we make the system be- 
havior look more and more classical. 



VI. SEPARABILITY CONSERVATION 

The last point to check is that our cloner does not 
create entanglement by itself, that is, it clones separable 
states into separable states. First, an important obser- 
vation is that our cloner is such that the input-to-single- 
clone transformation is a PPT map. Using Eqs. and 
(|ü|). and tracing the joint state \S)iz,a.,b,A over the ancilla 
A and one of the clones, say b, we arrive at the following 
expression for the state of the reference and the other 
clone 

Sn, a = \A\ 2 (t A <E> t B )n,a 

+ d 2 ((d 2 + 2)\A\ 2 + 2\B\ 2 + 4dVl(A*B)) 

x(\4>) a {4>\a <8> |0)b(^|b)tc,o 
+ d(d\B\ 2 + 2$t(A*B)) 

x(\(j>) A ((j>\ A (S>l B + t A (g)\(j)) B (<p\ B ) n , a (26) 

where TL A is the identity operator in the joint space of 
Alice's component of the reference 1Z and clone a, while 



\4>)a = d x/2 J2i=o \i) aMí) A,a is a ME state in this 
same space. (The same notations are used for Bob's ana- 
log quantities 1 B and \4>)b-) The cloning map is thus 
PPT since {Sh,o) Tb > 0, where T B stands for the par- 
tial transposition with respect to Bob's components of 
the reference 1Z and clone a. This PPT property ensures 
that the cloning of any isotropic state cannot increase 
its fidelity, hence its entanglement of formation ^l]]- In 
particular, all separable isotropic states are necessarily 
transformed into separable clones. 

In order to generalize this separability conservation 
property to all separable input states, outside the re- 
stricted class of isotropic states, we consider the cloning 
of a product state p A ® p B . By tracing (p A ® p B ) T Sn,a 
over the reference 1Z , we obtain for the first clone a state 
of the form 

Pa = \A\ 2 (1 A ® l B ) a 
+ ((d 2 + 2)\A\ 2 + 2\B\ 2 + Ad^R(A*B)) (p A <g> p B ) a 
+ (d\B\ 2 + 2^{A*B))( PA ®t B + t A ® PB ) a (27) 

where 1 A and 1 B are identities in Alice's and Bob's sub- 
spaces of clone a, respectively. Since all terms in l|27|l 
are product states and all coefficients are positive semi- 
defmitc for all d, we verify that p a is separable. By lincar- 
ity of the trace, this result also holds for any linear com- 
bination ^iPipf®pf withpi > and YliPi = 1)' ^hat is 
for a genèric separable state. Thus, we can conclude that 
our entanglement cloner transforms all initially separable 
states into separable clones. 

VII. CONCLUSION 

In conclusion, we have constructed an optimal (sym- 
metric) entanglement cloner, which is universal over the 
set of d x d ME states in arbitrary dimension d. On one 
hand, all separable input states are cloned by this cloner 
into separable states. In addition, the entanglement of 
isotropic states cannot be increased by the cloner (and 
we conjecture this property holds in general for any in- 
put state). On the other hand, the entanglement of the 
clones of ME input states is maximum. The optimiza- 
tion of the parameters of our QCM was performed by 
maximizing the fidelity of the clones, but the monotonic 
behavior of the entanglement of formation as a function 
of the fidelity for isotropic states guarantees that such an 
optimization maximizes the entanglement of the clones at 
the same time. We expect that entanglement is cloned 
"monotonically" for all (including non-isotropic) states, 
that is, higher entangled states result in higher entangled 
clones, and therefore the ME input states are those which 
generate the clones with the maximum achievable entan- 
glement. If this very natural assumption is right, then, 
based on our result, one can state that the maximal en- 
tanglement attainable by cloning is always below one half 
of the entanglement of the input state and saturates this 
value in the limit of large dimension d. This is consistent 
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with the idea that, since our QCM transforms separable 
states into separable clones, no additional cntanglement 
is produced by cloning, so we can only split the entan- 
glement of the input state between the two clones. This 
cxplains as well the asymptotic value of one fourth the 
initial cntanglement for the local cloner at the limit of 
large d. It is natural to expect all these conclusions to be 
vàlid for asymmetric entanglement cloners as well. 
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